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We simulate structural phase behavior of polymer-grafted colloidal particles by molecular Monte Carlo technique. 
Interparticle potential, which has a finite repulsive square-step outside a rigid core of the colloid, was previously con- 
firmed via numerical self-consistent field calculation. This model potential is purely repulsive. We simulate these model 
colloids in the canonical ensemble in 2 dimensions and find that these particles containing no interparticle attraction 
self-assemble and align in a string-like assembly, at low temperature and high density. This string-hke colloidal assem- 
bly is related to percolation phenomena. Analyzing the cluster size distribution and the average string length, we build 
phase diagrams and discover that the average string length diverges around the region where the melting transition line 
and the percolation transition fine cross. This result is similar to Ising spin systems, in which the percolation transition 
line and the order-disorder Une meet at a critical point. 

PACS numbers: 64.60.ah, 61.43.Er, 82.70.Dd 



I. INTRODUCTION 

Colloidal particles are known as hard particles made from 
glasses, metals, and other materials. Latex is an example of 
such artificial colloids. Natural colloids, e.g. pollen, are also 
seen in the environment. Colloids are known for their char- 
acteristic collective phenomena, for example glass transition 
and depletion interaction. Colloidal systems can be modelled 
by hard particle systems', owing to a sohd and undeformable 
characteristic of the particles. The phase behavior of hard par- 
ticle systems shows crystalization, which is called Alder tran- 
sition^'^. 

On the other hand, colloidal particles are liable to floccu- 
late and make deposition after a long time, which disturbs col- 
loidal dispersion. Polymers are, in industry, often grafted onto 
the colloidal particles to preclude the deposition and stabi- 
lize the dispersion. These are called polymer-grafted colloidal 
particles. We can control the interaction between polymer- 
grafted colloids by adjusting the length, species, the grafting 
density, and the chemical quality of the grafted polymers'*"**. 
Such usefulness of polymer-grafted colloids brings various 
applications and plays an important role in industry. 

Simulating a simple model system via particle Monte Carlo 
simulation technique, we study the efi'ect of the grafted poly- 
mers on the phase behavior of colloidal systems. In the 
present article, we show that particles interacting via spher- 
ically symmetrical repulsive square-step pair potential, which 
contains no attraction, self-assemble and finally align in 
strings. It is also shown that this string-like assembly indicates 
a similarity to percolation transition and critical phenomena^. 

A pair of polymer-grafted colloids interact with, in addition 
to their rigid cores, the grafted polymers outside these rigid 
cores. For example, let us consider an aqueous solution of the 
colloidal particles, diblock copolymers consisting of a long 
hydrophobic block and a short hydrophilic block are grafted. 
The hydrophobic blocks of these grafted diblock copolymers 
form a melt polymer brush, while the hydrophilic blocks form 
a thin shell swollen by the solvent. The rigid core is covered 



in the melt brush, which is covered in the swollen shell. 

A pair of swollen brushes repels each other owing to the 
excluded volume effect'". On the other hand, slight attractive 
interaction potential is found between a parr of melt brushes 
upon contact'", since the polymer chain conformation relax 
when the two melt polymer brushes contact. The so-called 
Derjaguin approximation yields the effect of the curvature of 
the grafting surfaces on the interaction between grafted poly- 
mer brushes", provided that the radius of the curvature is far 
larger than the brush thickness. However, as long as the ra- 
dius of the curvature is the same order of the brush thickness, 
no good analytical approximation is known. Numerical self- 
consistent field (SCF) calculation^^''" is necessary for such a 
situation. 

Our SCF calculation has shown that, with some parameter 
sets, the interaction between a pair of our polymer-grafted col- 
loidal particles reduces to a repulsive step potential containing 
no attraction'. We approximate this potential by spherically 
synnmetrical square-step potential with a rigid core, 

0(r) = oo r < <T\, 

0(r) = eo (> 0) <Ti < r < ctj, 
4>(r) = o"2 < r, 

where r denotes the distance between the centers of the pair 
of the particles, cri denotes the diameter of the colloids, and 
0-2 is the diameter of the outer core formed by the polymer 
brush. Note that the height of the repulsive potential step eo, 
which stems from the conformational entropy of the grafted 
polymers, is dependent on the temperature'. 

Due to the repulsive square-step in 0(r), the phase behavior 
of the present system depends on temperature. This is differ- 
ent from the behavior of the hard particle systems^'^. Boltz- 
mann factor exp(-eo/^BT) determines the probability that a 
pair of particles get on the step in the potential, where de- 
notes the thermal energy. The phase behavior at extremely low 
or high temperature is easily understood, although our model 
potential, is invalidated in such extreme regions of the 
temperature. 
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At extremely high temperature, the phase behavior of our 
system is identical to the behavior of a system of hard spher- 
ical particles with diameter cri since the height of the step, 
eo, is vanishingly small compared with the thermal, or kinetic, 
energy'. 

At extremely low temperature, the repulsive step of (p{r) 
becomes far higher than the thermal energy. Due to this ex- 
tremely high potential energy barrier, phase diagrams of the 
system is consistent with those of a system of hard spheri- 
cal particles with diameter 0-2, when the system volume, V, 
is larger than the close-packed volume of the outer core, Vq- 
The pressure suddenly rises when V is reduced across Vo, i.e. 
the compression of the system box. This results in a phase 
transition at V = Vq, which is similar to the Alder transition'. 

In early works, a step potential with a rigid core or one with 
both a rigid core and a square-well were introduced as a model 
potential for the purpose of qualitatively simulating anomalies 
of phase behavior, e.g. anomalous behavior of the melting 
curves and the critical behavior of phase diagrams, observed 
in systems of water, cesium, cerium, and others'^"'*. 

At finite temperature T, the model system forms dimers, 
lamellae, and other various assembling structures'"'^"^^. Glass 
transition^-' is observed in the same model system. These as- 
semblies are also found in a system of particles interacting via 
continuous repulsive potential similar to our model potential 
<p(r)^'^. A variety of ground states of the same model system 
have been discovered at zero temperature via genetic algo- 
rithms^^. In the present study, we simulate our model system 
for the sake of constructing complete phase diagrams at finite 
T ml dimensions. Simulation results in 3 dimensions will be 
published in our forthcoming article^^. 



II. SIMULATION METHODS 

Monte Carlo simulations are performed in 2 dimensions 
in the canonical ensemble (A^VT-ensemble) via the standard 
Metropolis algorithm^^-'^^. and V denote the number of par- 
ticles and the system volume (area) respectively. Mersenne 
Twister is chosen as a random number generator for our sim- 
ulations^'"^ ^ In the initial state the particles are placed on 
triangular lattices in the system box with periodic boundary 
conditions. In one simulation step, a particle is picked at ran- 
dom and given a uniform random trial displacement within a 
square whose sides have length Q.4cr[. A Monte Carlo step 
(MCS) is defined as trial moves, during which each parti- 
cle is chosen for the trial displacement once on average. After 
5.0x 10^ MCS, by which the system relaxes to the equilibrium 
state except at low temperature and high density, we acquire 
data every 10^ MCS and gel 100 independent samples of parti- 
cle configurations. Simulation results at A' = 1200 are mostly 
shown in the present article. We have confirmed, however, 
that the physical properties of the simulation system does not 
significantly change even when we run simulation for a larger 
system with N = 4800. 

(Ti and eo are taken as the unit length and the unit energy 
respectively. Vq (Vq - V3/2 in 2 dimensions) is used as 
the unit of the volume (area). Dimensionless temperature is 



defined asksT/eo. 



III. SIMULATION RESULTS AT o-j/o-i = 2.0 

Here we discuss simulation results at 0-2/0-1 = 2.0, a typical 
case for the strongly curved grafted-brush. 

As a particular example, we show a snapshot of the system 
at = 10800, V/Vq = 0.7, ksT/eo = 0.1, and 5.5 x 10*' MCS 
in Fig. 1. We find that the colloidal particles self-assemble 
and finally align in strings. We called this efl^ect "frogspawn 
effect" owing to the characteristic particle configuration and 
called this structure "string-like assembly"' ""^'. This eff'ect 
has experimentally been confirmed lately^^. This string-like 
assembly in 2 dimensions is also observed at 0-2/0-1 - 2.0^^ 
or 2.5'^, where the system is cooled from high T to low T, and 
in a model system that consists of particles interacting via con- 
tinuous repulsive potential similar to (f>(r) in 2 dimensions^'*. 
We have discovered that the string-like assembly also appears 
in our model system in 3 dimensions', which will separately 
be discussed in our forthcoming paper^*. 

A. Definitions of string length and cluster size 

For the purpose of characterizing the structure of the parti- 
cle aggregation, string length and cluster size are defined. 

First, a "bond" is defined as an overlap between a pair of 
particles. A single network of bonds consisting of M particles 
is defined as a "cluster with the size M". A particle that has 
more than two bonds in a cluster is defined as a "junction" of 
the network. 

As is shown in Fig. 1, there are clusters with the same 
size and the different configurations of the network. We in- 
troduce "string length" for quantitative distinction between 
the network configurations of clusters. A string is defined 
as a sequence of particles, each of which has exactly two 
bonds. The string terminates with junctions or particles with 
only one bond unless a circular configuration is composed. 
String length is defined as the number of particles along the 
string, including both the ends. A circular cluster consisting 
of three particles are excluded from the definition of the string 
although all the particles in such a cluster have two bonds. 



B. Stability of tlie string-like assembly 

Here we check the dependence of the final colloidal as- 
sembly on the initial assembly. For this purpose, simulations 
with the same set of parameters as Fig. 1, i.e. V/Vq = 0.7 
and ksT/sQ = 0.1, have been performed starting from 4 dif- 
ferent homogeneous and/or inhomogeneous initial configura- 
tions, which are perfect string-Uke assembly aligned in one 
direction, the close-packed configuration of the inner cores, 
and two other configurations. The string length distributions 
and the cluster size distributions, which characterize the par- 
ticle configuration found in the system, are analyzed for these 
4 simulation runs with the same parameter set. The analyzed 
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distributions indicate that all the simulation runs starting from 
different initial assembly result in the same string-like struc- 
ture. This result demonstrates the stability of the string-like 
assembly. On the other hand, these cluster size distributions 
show an exponential distribution with the cluster size, which 
is similar to the formation of clusters in percolation systems 
in its non-critical region. This correspondence is discussed in 
section HID. 

Considering the effect of the internal energy, we investi- 
gate the reason for the string-like assembly. At extremely low 
T, the particles tend to get down the step of the potential to 
decrease the internal energy. In regions of V/Vo < 1.0, how- 
ever, enough free volume to allow all the particles to get down 
the step is not found in the system. Therefore, at extremely 
low T and V/Vo < 1.0, the system builds long strings since a 
string-like linear cluster takes lower internal energy than a ho- 
mogeneous crystalline structure, which is chosen as the initial 
configuration of our simulation. As the temperature increases, 
the string length decreases due to large probablity of overlaps 
between the particles at high T . Long strings are absent at 
y/Vo ~ 1.0 since the particles seldom overlap. At high den- 
sity, since the small free space rarely allows the particles to 
get down the step, the particles build only a few strings and 
most particles are junctions. 



C. Average string length at crijai = 2.0 

For the sake of locating the region of the string-like assem- 
bly, we analyze the average string length which is presented in 
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FIG. 2. The average string length at (T2lo'i = 2.0 and A' = 1200. 
Blue regions represent short strings and red ones long strings. Grids 
denote the points where the simulation is performed. 



Fig. 2. The average string length has a sharp peak in a region 
of 0.1 < ksTjeo < 0.2 and 0.6 < V/Vo < 0.7. Long strings 
are observed only in the vicinity of this peak. At high density, 
there are only a few strings since particles are rarely allowed 
to get down the potential step from the initial crystalline con- 
figuration. At y/Vo > 1-0 and low k^T jeo, particles seldom 
overlap and most particles are isolated. 
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FIG. 3. Occurrence probability of percolated clusters for the same 
system as Fig. 2. (rijai = 2.0 and A' = 1200. 




FIG. 4. Cluster size distribution at crj/cri = 2.0, VjVo = 0.98, 
kgT jeo = 1.0, and N = 4800. This parameter set is located on the 
percolation transition line in Fig. 3. log(i)-log(«(.s)) graph is plotted. 
A solid black line shows a linear fitting of this graph, y = -r * x + b 
withT = 1.9 andZ? = -2.3. 



The graph of the average string length with a resolution 
finer than Fig. 2 indicates that the maximum average string 
length exceeds 100 at V/Vo = 0.64 and kBT/eo = 0.12, which 
suggests a divergence of the string length''. Simulation results 
via 3-reservoirs Monte Carlo method, which show parallel- 
arranged globally-anisotropic defect-free string-like assem- 
bly, corroborate this divergence^^'^"*. 3-reservoirs method al- 
lows us to efficiently simulate the equilibrium state. This re- 
sult via 3-reservoirs method is discussed in section HIE. 



D. Percolation phenomena at 0-2/0-1 = 2.0 

As is observed in Fig. 3, at = 1200 and kBT/eo = 0.1, the 
maximum cluster size suddenly changes between V/Vo = 0.6 
and V/Vo = 0.7. At V/Vq = 0.7 and 0.8 the maximum cluster 
size found in the system equals no more than approximately 
70 in our simulation runs, whereas at V/Vo = 0.6 and 0.5 ex- 
ceeds 1000. At V/V{) > 0.7 many small clusters are observed 
in the system, whereas at V/Vo ^ 0.6 most particles in the 
system compose one large cluster In order to understand this 
sudden change, we calculate the probability of the occurrence 
of a large cluster that bridges both the sides of the system 
box, i.e. a percolated cluster^^. The occurence probability of 
the percolated cluster at each hgT/eo and V/Vo is presented in 
Fig. 3. A sharp boundary around V/Vo ~ 0.7-1.0 divides per- 
colation and non-percolation regions. The occuiTence proba- 
bility of the percolated clusters changes abruptly at this sharp 
boundary, which suggests a similarity of this phenomenon to 
the usual percolation phenomena^^. 

At low density, V/Vq > 1.0, the particles do not frequently 
overlap due to enough free volume to allow all particles to get 
down the interaction potential step. This results in the bound- 
ary in Fig. 3 asymptotically approaching to a line V/Vo = 1.0 
at high kgT/eo- On the other hand, at low hgT/eo, the bound- 
ary is located at V/Vo ~ 0.65, between Alder transitions of 
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FIG. 5. Cluster size distribution at iTt/o-] = 2.0, V/Vq = 0.7, 
kgT /ef) = 0.1, and N = 4800. This parameter set is located outside 
the percolation transition line in Fig. 3. i-log(ii(i')) graph is plotted. 
A solid black line shows a linear fitting of this graph, y = -a * x + b 
with a = 0.16 and b = -3.7. 



outer cores and inner cores. 

Here we study the relationship between our string-like as- 
sembly and the percolation phenomena. We calculate the 
cluster size distribution, n{s), along the boundary curve in 
Fig. 3. Here, the cluster size distribution is defined as, n{s) = 
M{s)/N, where s denotes the cluster size and M{s) is the num- 
ber of clusters with the size s found in the system. An example 
of n{s) on the percolation transition line in Fig. 3 is presented 
in Fig. 4 in double logarithmic scales, which shows a relation, 
n{s) oc 5 with T - 1.9, where t means the Fisher exponent, a 
critical exponent, of the percolation transition of our system. 
This power law with t - 1 .9 is unchanged along the percola- 
tion transition line. Outside the percolation transition line, an 
exponential distribution is observed, which is also similar to 
the usual percolation phenomena. An example of this expo- 
nential distribution is plotted in Fig. 5. 

Divergence of the string length in Fig. 2, an anomalous 
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FIG. 6. The probability density of the square diaplacement of the 
particles, Wsd, during 1 x 10^ MCS at cr2la-i = 2.0, A' = 1200, and 
kgT lea = 0.01. The result at V/Vo = 1-2 is plotted by a grey solid 
line and V/Vo = 1.3 by a black broken line. The number of analyzed 
particle configurations, i.e. the number of samples, is equal to 9. 



phase behavior, is found around the percolation transition line. 
This result indicates a similarity between our string-like as- 
sembly and critical phenomena. One of the critical exponents, 
T, has been determined above in the present section, whereas, 
owing to the insufficient resolution of the data, other critical 
exponents are not evaluated in the present work. 



E. Thermodynamic phase diagram at 0-2/0"! = 2.0 

A thermodynamic phase diagram is constructed in the 
present section. 

As has been discussed in section I, at extremely low tem- 
perature and y/Vo > 1.0, the phase behavior of our system 
coincides with the behavior of the hard particle system of the 
outer cores. Alder transition of the hard particle system occurs 
at y/Vo ~ 1.3 in 2-dimensions^. For the purpose of confirm- 
ing this Alder transition of the outer cores of our colloids, in 
low temperature regions ksT jeQ < 0.1, the probability den- 
sity of the square displacement of the particles is analyzed at 
W^o = 1-2 and 1.3. The results, e.g. presented in Fig. 6, 
significantly change between V/Vq -1.2 and 1.3, which in- 
dicates the Alder transition. In other regions of V/Vo and 
ksT/eQ, a thermodynamic phase diagram is constructed based 
on the mean-square displacement of the particles (MSD) dur- 
ing 99 X 10^ MCS. MSD at V/Vq = 1.3 and /ter/eo = 0.01 is 
chosen as the standard point, owing to the Alder transition of 
the outer cores. When MSD is smaller than this standard, the 
system is regarded as a solid. When MSD is larger than the 
standard value, the system is in a fluid phase or a coexisting 
phase of both the phases. 

The constructed thermodynamic phase diagram is given in 
Fig. 7. At high ksT/eo Alder transition of the inner cores 
of the colloids occurs in regions of 0.3 < V/Vq < 0.4, 
whereas Alder transition of the outer cores, i.e. the transition 
at V/Vq * 1.3, does not. At low ksT/eo only Alder transi- 
tion of outer cores occurs. At intermediate ksT/eQ, the solid 
melts when not only V/Vq but also kgT/ea increases. These 
results are consistent with the phase behavior discussed in the 
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FIG. 7. Thermodynamic phase diagram at 0-2/0-1 = 2.0 and N = 
1200. Dark grey regions represent solid phases and light grey fluid 
or coexistence of both the phases. Alder transition points of the inner 
and the outer cores are pointed by arrows above the panel. 



introduction, section I. 

The region where the string-like assembly, as is shown in 
Fig. 1, and the percolation transition are located, i.e. the re- 
gion of 0.1 < keT/eu < 0.2 and 0.6 < V/Vo < 0.7, is found 
in the solid phase between the two Alder transitions. This 
means that the string-like assembly is an amorphous solid, a 
non-equilibrium state. 

Equilibration in the canonical ensemble demands a long 
compuational time. Therefore, the equilibrium state has re- 
cently been simulated via 3-reservoirs method, which pro- 
vides shortcuts in the phase space from non-equilibrium to 
equilibrium states-^-'--''*. Moreover, 3-reservoirs method allows 
the system itself to finely and simultaneously tune and the 
system box size to the equilibrated structure. Unlike other 
simulation methods, due to this unique feature, anisotropy of 
the equilibrated ordered structure is kept in a global scale. The 
periodicity of the structure is not constrained to the fixed sys- 
tem box size. Defects are removed from the ordered struc- 
ture. The simulation results via this 3-reservoirs method have 
shown the globally-anisotropic defect-free string-like assem- 
bly as the equilibrated structure-'^'^"*, presented in Fig. 8. 



IV. SIMULATION RESULTS AT VARIOUS STEP WIDTHS 

We have found the string-like assembly at 0-2/a-i = 2.0. 
Here we study a range of (T2/0-1, where the string-like assem- 
bly is observed. 

When the step width 0-2/0-1 is taken to a limit, 0-2/0-] — > 1, 
the potential step disappears and the phase behavior of the sys- 
tem is equivalent to the behavior of the hard particle system 
with a particle diameter crj . This means that the string-like 
assembly is absent in this limit. On the other hand, in a limit 
of 0-2/0-1 — » 00, the inner hard core vanishes, which also re- 
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FIG. 8. Snapshot of the system simulated via 3-reservoirs method^^'. 
Black dots represent the centers of the particles and grey lines denote 
the bonds between the particles. hgT leo = 0.12, Pa-j/eo = 1.1, and 
yj' = 29.93. P and fi' denote the pressure and the dimensionless 
chemical potential respectively. 

suits in the absence of the string-like assembly. Therefore, 
the string-like assembly is observed only in an intermediate 
region of step width 0-2! cr\. This result contradicts a conjec- 
ture by Malescio and Pellicane'** that this assembly should be 
found at any step width satisfying a condition 0-2/0-1 ^ 1. 

We simulate the system at various step widths and deter- 
mine this intermediate region of cr2lo-\ in the present section. 
The simulation results in section III indicate this intermediate 
region ranges around 0-2/0-1 - 2.0, which simulation results 
in the present section corroborate. Other physical properties 
of the system at various step widths are also studied. 

A. Average string length at various step widths 

The average string length at 0-2/0-1 - 1.1, 1.5, 1.9, 2.5, and 
3.0 are given in Fig. 9. At cr2lo-\ - 1.9 and 2.5, the string- 
like assembly as is shown in Fig. 1 is observed, whereas, at 
0-2/0-1 = 1.1, 1.5, and 3.0, the long strings are not observed 
in our simulation. The same graph as Fig. 9(b), but with high 
resolution, is presented in Fig. 10, which demonstrates that the 
resolution of the graphs does not affect this result and that the 
divergence of the string length is absent at this 0-2/0-1. As an 
example, a snapshot of the system on a peak of Fig. 10 is given 
in Fig. 11. Our simulation results indicate the region of the 
string-like assembly ranging in an interval 1.7 % 0-2I0-1 g 2.8, 
i.e. in the vicinity of cr^lcrx - 2.0. 

The boundary of the step width, cr^l^'x - 2.0, divides 
characteristics of the local particle configuration. Below this 
boundary, i.e. in the region of 0-2/0-1 < 2.0, the phase behav- 
ior of the system is similar to the phase behavior at the limit 
of 0-2/0-1 — » 1. When the outer cores of a pair of the particles 
are in close contact as is sketched in Fig. 12(a), the 3rd par- 



ticle is disallowed to be placed on a straight line between the 
pair, due to overlaps of the inner cores. 

On the other hand, at this boundary of the step width 
0-2/0-1 - 2.0, the 3rd particle is allowed to be arranged on 
this straight line, as is illustrated in Fig. 12(b). In addition, 
in this particle configuration, the inner cores of the 3 particles 
are in close contact. The displacement of the 3rd particle from 
this position is limited to the direction normal to this straight 
line. The 3rd particle is pinched and confined by the other 
2 particles. When particles are added to these 3 particles and 
closely packed along the straight line (Fig. 12(c)), not only the 
inner cores but also the outer cores of all the particles, exclud- 
ing both the end particles, simultaneously pinch and confine 
each other in the same way. This means that the inner core 
is pinched by the nearest neighbors and that the outer core by 
the next nearest neighbors. Therefore, the string-like assem- 
bly is stable in the direction of the straight line. When these 
straight strings are closely arranged in the same direction, the 
string-like assembly is also stable in the lateral direction since 
an additional energy cost is required to move and force one 
particle of a string into the adjacent string. Therefore, the 
system tends to keep the linear, i.e. string-like, assembly at 
0-2/0-1 = 2.0. 

For the case, 0-2/0-1 > 2.0, the 3rd particle is able to be 
arranged on the straight line between the pair of the particles 
(Fig. 12(d)). However, due to the free space between the pair, 
the 3rd particle is allowed to freely move. This disrupts the 
string-like assembly. When the inner cores of the 3 particles 
are closely packed along the straight line, this particle config- 
uration demands eg additional energy cost, which also disrupts 
the string-like assembly. 

In 3 dimensions, additional degrees of freedom reduce the 
pinch and the confinement of the particles and decrease the in- 
terval of 0-2/0-1 where the string-like assembly is observed^^. 
In other words, in 3 dimensions, the string-like assembly is 
observed only in the close vicinity of 0-2/0-1 - 2.0. 



B. Percolation phenomena and thermodynamic phase 
diagrams at various step widths 

The occurrence probability of percolated clusters and ther- 
modynamic phase diagrams at 0-2/0-1 - 1.1, 1.5, 1.9, 2.5, and 
3.0 are presented in Figs. 13 and 14. These are qualitatively 
consistent with the results at 0-2/0-1 = 2.0. The Fisher expo- 
nent, T - 1.9, is also kept at these step widths. These results 
indicate that the discussion on the results at 0-2/0-1 = 2.0 also 
apply to the results at these 0-2/0-1. 

In addition to the simulation results at 0-2/0-1 = 2.0, the 
simulation results at 0-2/0-1 - 1.9 and 2.5 show that the aver- 
age string length diverges around the region where the melt- 
ing transition line and the percolation transition line cross^. A 
similar behavior is found in Ising spin system where the per- 
colation transition line and the order-disorder line meet at the 
critical point-'^. 
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FIG. 9. The average string length at = 1200. ctt/o-i = (a): 1.1, (b): 1.5, (c): 1.9, (d): 2.5, and (e): 3.0. Red regions represent long strings 
and blue short strings. The average string length at (Tilai = 2.0 is given in Fig. 2. 
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FIG. 10. The average string length at A' = 1200 and 0-2/0-1 = 1.5. 
The same graph as Fig. 9(b), but with high resolution. The long 
strings and the divergence of the string length are absent at this 
0-2/0-1. 




FIG. 11. Snapshot of the system on a peak of Fig. 10, i.e. at 0-2/0-1 = 
1.5, N = 1200, V/Vo = 0.84, ytgr/eo = 0.06, and 5.0 x 10* MCS. 
Black dots represent the centers of the particles and grey lines denote 
the bonds between the particles. Long strings, as is shown in Fig. 1, 
are not observed. 



V. CONCLUSIONS 

Utilizing Monte Carlo simulation technique, we have stud- 
ied the phase behavior of 2-dimensional systems of hard par- 
ticles with square-step repulsive interacting core. The string- 





FIG. 12. Sketch of local configurations of a few particles arranged 
on a straight line, (a): 0-2/0-1 < 2.0. (b) and (c): 0-2/0-1 = 2.0. (d): 
0-2/0-1 > 2.0. 



like assembly of particles between the two Alder transition 
densities, each of which corresponds to the inner hard core 
and the outer hard core, has been found. This string-like 
assembly is related to the percolation and the critical phe- 
nomena. Our model system has recently been studied in ex- 
periments, in which the string-like assembly has been con- 
firmed^^. 

By simulating the system at various step widths, 0-2/0-1, we 
determined dependence of these results on the step width. We 
have discovered that the string-like assembly, as is shown in 
Fig. 1, is observed in the interval 1.7 g 0-2/0-1 g 2.8, i.e. in 
the vicinity of 0-2/0-1 = 2.0. At this step width 0-2/0-^ = 2.0, 
when the particles are closely and linearly packed on a straight 
line, both the inner cores and the outer cores of the particles 
simultaneously pinch and confine each other. Therefore, the 
region of the stiing-like assembly ranges around 0-2/0-1 - 2.0. 

The Fisher exponent of the system is found to be t = 1.9 at 
any step width 0-2/ o-\. 

Real systems coiTesponding to the model particles, which 
are interacting via the square-step repulsive potential, (/)(/), 
or the square-step potential with the attractive square-well, 
were absent in previous studies'^"'^-'^-'^. However, in the 
present study, the interaction potential (/>(r) has been linked 
to polymer-grafted colloidal particles, i.e. real particles, and 
been derived from the numerical self-consistent field (SCF) 
calculation^. 

With the use of the SCF method, we are able to calculate 
the interaction potential between a pair of polymer-grafted 
colloidal particles for various physical and chemical param- 
eter sets, e.g. polymer chain length, block-ratio, grafting den- 
sity, monomer species, and solvent quality. These results are 
fed into the particle Monte Carlo scheme. Thus, combin- 
ing our particle Monte Carlo and SCF calculation results, we 
can predict and control the phase behaviour of any diblock 
copolymer-grafted colloidal systems. This methodology al- 
lows us to extend this approach to more sophisticated block 
copolymers and colloidal systems. 

The string-like assembly itself is also useful for application. 
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FIG. 13. Occurrence probability of percolated clusters a.tN= 1200. <T2/(T| = (a): 1.1, (b): 1.5, (c): 1.9, (d): 2.5, and (e): 3.0. Corresponding 
occurrence probability of percolated clusters at (T2lcri = 2.0 is given in Fig. 3. 
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FIG. 14. Thermodynamic phase diagrams at A' = 1200. crilo'i = (a): 1.1, (b): 1.5, (c): 1.9, (d): 2.5, and (e): 3.0. Corresponding 
thermodynamic phase diagram at ailcri = 2.0 is given in Fig. 7. Alder transition points of the inner and the outer cores are pointed by arrows 
above the panel. 
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As an example, when the hydrophilic blocks of the grafted 
polymers that cover the colloid with a thin swollen shell are 
composed of dielectric material and the thick hydrophobic 
brush that fills the inside of the shell consists of conductive 
material, our system can serve as a nano- switching device; 
the system can conduct the electricity only in the percolated 
phase. Another example is display devices with use of the dif- 
ference between the string-like assembly and the disordered 
phase. 

The string-like assembly is also found in 3-dimensional 
systems'. A fuU report will be published in our forthcoming 
paper^*. 
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